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Abstract 
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The most general admissible boundary conditions are derived for 
an idealised Aharonov-Bohm flux intersecting the plane at the origin 
on the background of a homogeneous magnetic field. A standard tech- 
nique based on self-adjoint extensions yields a four-parameter family 
of boundary conditions; other two parameters of the model are the 
Aharonov-Bohm flux and the homogeneous magnetic field. The gen- 
eralised boundary conditions may be regarded as a combination of the 
Aharonov-Bohm effect with a point interaction. Spectral properties of 
the derived Hamiltonians are studied in detail. 

1 Introduction 

The purpose of this paper is to determine the most general admissible bound- 
ary conditions for the Aharonov-Bohm (AB) effect in the plane on the back- 
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ground of a homogeneous magnetic field, and also to investigate the basic 
properties of Hamiltonians obtained this way. The history of the effect goes 
four decades back and starts from the observation of Aharonov and Bohm 

that the behavior of a charged quantum particle is influenced by a mag- 
netic flux even if the field is zero in the region where the particle is localized. 
A particularly elegant treatment is possible in case of an idealized setup in 
which the AB flux is concentrated along a line perpendicularly intersecting 
the plane, conventionally at the origin 

The boundary conditions of the last mentioned paper are not the most 
general ones; the full family of such conditions giving the AB effect in the 
plane was derived in fl7|, and simultaneously also in ||. These generalised 
boundary conditions may be interpreted as a combination of the AB effect 
with a point interaction supported, too, at the origin, although this is just 
one possible point of view. In any case they can be described and investigated 
by the technique of self-adjoint extensions which is in principle the same one 
as that used in the paper [|] in which two-dimensional point interactions 
were introduced. 

A natural question is what happens if such a system is placed into a back- 
ground homogeneous magnetic field. This problem attracted some attention 
recently, even with a controversy: the papers ||, |, || consider the "pure" AB 
effect in this setting for the Pauli operator, i.e. a spin 1/2 particle. The last 
named property leads to specific behavior related to the Aharonov-Casher 
effect, which we will not discuss here. 

Our aim here is different: we are going to consider a spinless particle 
with a point flux and a homogeneous background, and ask about the most 
general class of boundary conditions analogous to those of |7], |J. The basic 
difference between the situations without and with a homogeneous magnetic 
field is that in the former case the spectrum is absolutely continuous and 
equal to the positive half-line possibly augmented with at most two negative 
eigenvalues (depending on the choice of boundary conditions) while in the 
latter case the spectrum is pure point and the point flux and interaction gives 
rise to eigenvalues in each gap between neighboring Landau levels. Our goal 
is to discuss these spectral properties in detail. 
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2 Formulation of the problem, preliminaries 



We consider the symmetric operator 

L = -(V-A(V)) 2 , Dom(L) = C -(R 2 \{0}), 

where the vector potential A is a sum of two parts, A = A^f + A^, with the 
part A hmf corresponding to the homogeneous magnetic field in the circular 

gauge, 

iB 

^hmf = ^" (-^2^1 + Xldx 2 ), 

and with the part Aab corresponding to the idealised AB effect, 

Aab = w^(-x 2 dx 1 + x 1 dx 2 ), r 2 = x 2 + x 2 2 . 
27rr z 

Without loss of generality we may assume that B > 0. Further, we rescale 
the Aharonov-Bohm flux, 

to have a variable which expresses the number of flux quanta and, as usual, we 
make use of the gauge symmetry allowing us to assume that a G ]0, 1[ . Hence 
the case $ G 2nZ is excluded since it is gauge equivalent to the vanishing 
AB flux. Our goal is to describe all the self-adjoint extensions of L as well 
as to investigate their basic properties. 

It is straightforward to determine the adjoint operator L* , 

■ip G Dom(L*) ^GL 2 (R 2 ,d 2 x)niJ 1 2 o ' c 2 (M 2 \{0}) 

and (V - A(V)fi) G L 2 (R 2 , d 2 x). 

Next we can employ the rotational symmetry when using the polar coordi- 
nates (r, 6) and decomposing the Hilbert space into the orthogonal sum of 
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the eigenspaces of the angular momentum, 

L 2 (R 2 ,d 2 x) = ^®L 2 (R + ,rdr)®Ce im9 . 



In the polar coordinates the operator L (and correspondingly L*) takes the 
form 

1 „ „ 1 / „ Br 2 ^ 2 



L — — d r rd r H — - — ide + a + 



r r 2 \ 2 

The operator L* commutes on Dom(L*) with the projectors P m onto the 
eigenspaces of the angular momentum, 

P m ^(r, 0) = J- ^(r, 0') d6' , 

^ Jo 

and therefore L* decomposes in correspondence with the orthogonal sum ([[]), 

L * = E Vw (2) 

mgZ 

Thus we can reduce the problem and work in the sectors RanP m , mGZ. 
For a given spectral parameter A G C we choose two independent solutions 
(except of particular values of A) of the differential equation 

1 1 / Br 2 \ 2 \ 

-d r rd r + — {m + a + ^-j I g(r) = Xg(r), (3) 



namely 

^(A;r)=rl m+a lW/3(m,A), 7 (m), 



2 

9m 



(A;r) =r |m+a| G (/5(m,A),7(m) 



Br 2\ 


) exp ( 
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where 

(3(m, A) = - ( 1 + m + a + \m + a\ — — 
7(m) = 1 + \m + a\ . 
Here F and G are confluent hypergeometric functions ||, Chp. 13], 



(5) 



F((3,j,z) = J2 



n=0 n! 



and 



G ^ 7, *) = rff- 7 +i) F(A 7 ' 2) + zl ' 1F{(3 ~ 7 + x ' 2 " 7 ' 2) ' 

(6) 

Notice that 7, 2) and 7, z) are linearly dependent if and only 
if /? G — Z+. Moreover, F(/3,j,z) is an entire function, particularly, it is 
regular at the origin while G(/3, 7, z) has a singularity there provided 7 > 1 
and /3 ^ — Z + , and in that case it holds true that 

lim z^Gtf^z) ih-i) 



Thus in the case when 1 < 7 < 2 we have the asymptotic behaviour, as 

We shall also need some information about the asymptotic behaviour at 
infinity. When z — > +00 it holds true that 



m 7, *) = r^Tj) (-*)-' (1 + (^)) + e ^' 7 (1 + ) 



(8) 
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and 



G(f3, 1 ,z) = z-P(l + 0(z- 1 )). 

3 The standard Aharonov-Bohm Hamiltonian 

With the above preliminaries it is straightforward to solve the spectral prob- 
lem for the standard AB Hamiltonian as we mentioned in the introduction. 
This means to solve the eigenvalue problem 

L*ip = Xip 

with the boundary condition 

lim if>(r, 6) = 0. (9) 

r— »0+ 

By virtue of the decomposition (||) the problem is reduced to the countable 
set of equations 

(L*) m f = Xf,me Z, 

and hence to the differential equations (|3|). 

The solution <7^(A; r) of (|3|) is ruled out because it contradicts the con- 
dition (|) and the solution g>^(A;r) belongs to L 2 (IR + ,rdr) if and only if 
P(m, A) = — n, with n G Z + . Since it holds 

F(-n, l + a,z) = ^'^^W K{z), n G Z+, 
T(n + a + 1) 

we get a countable set of eigenvalues, 

A m , n = B (m + a + |m + a\ + 2n + 1), m G Z, n G Z + , 
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with the corresponding eigenfunctions 




( 



Br 2 



Br 2 




n 



2 



exp 
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where 




±(|m+a|+l) 



1/2 



m,n 



7T F(n + \m + a\ + 1) 



are the normalisation constants. 

As it is well known if we fix m E Z then the functions {f m ,n( r , $)}^o 
form an orthonormal basis in L 2 (IR + , r dr) ®C e im0 and so the complete set of 
eigenfunctions {/ m ,„(r, 6 l )} mg ^, «ez+ is an orthonormal basis in L 2 (R + , r dr)® 
L 2 ([0, 27r], d8). Since all the eigenvalues \ m ,n are real we get this way a well 
defined self-adjoint operator which is an extension of L. We convention- 
ally call it the standard AB Hamiltonian and denote it by H AB . Thus the 
spectrum of H AB is pure point and can be written as a union of two parts, 



a(H AB ) = a pp (H AB ) = {B(2k + 1); k E Z + } U {B(2a + 2k+l); k E Z + }. 



Notice that the eigenvalues belonging to the first part are nothing but the 
Landau levels. All the eigenvalues B(2k + 1) have infinite multiplicities while 
the multiplicity of the eigenvalue B {2a + 2k + 1) is finite and equals k + 1. 

A final short remark concerning the Hamiltonian H AB is devoted to the 
Green function. Naturally, the Green function is expressible as an infinite 
series 



G AB (z;r l ,6 1 ,r 2 ,9 2 ) = ±- £ G^(z; n, r 2 ) e™<*-*> 



m=— oo 
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where 

> v |m+a|+l 



Gi B (z;n,r 2 ) = 2(-j (r^)^ exp B{r* + r 2 2 ) J 

oo 

E 



oo 

71 

X 



n T(n + \m + a\ + 1) 

n=0 v 1 1 ' 



i?(m + a + \m + a\ + 2n + 1) — z 



The radial parts can be rewritten with the aid of the standard construction 
of the Green function for ordinary differential operators of second order, 

/ fi\\m+a\+l / 1 

Gi B (z;n,r 2 ) = l^-j (r 1 r 2 )l m+ "lexp(^- i 5(r 1 2 + r 2 2 ) 

r( — w(m, z)) 

x wi — , — rrrr F (- w ( m , z),\m + a\ + 1, r<) 

r(|m + a| + 1) 
x G{— w(m, z), \m + a\ + 1, r>) 

where 

z 1 

w(m, z) = — — (m + a + |m + a\ + 1) 

and r< = min(ri,r 2 ), r> = max(ri,r 2 ). This amounts to the identity 

n! LJfeO L" n (y 2 ) 



E 



T(n + (7 + 1) n — w 
1 (<7 + 1J 

We do not expect that a simpler form for the Green function could be derived 
since the Hamiltonian H AB enjoys only the rotational symmetry. 
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4 Self-adjoint extensions of L 



Recalling what has been summarised in Section ^| it is easy to determine the 
deficiency indices. The solution ^(±?;r) diverges exponentially at infinity 
(cf. (H)) while g^ n (±t;r) behaves well at infinity but has a singularity at 
the origin of the order r -\ m+a \. Thus gr^(±i; r) G L 2 (M + ,rdr) if and only if 
m = — 1 or m = 0. This means that the deficiency indices are (2,2). For a 
basis in the deficiency subspaces J\f± % we can choose 

{/m,±(r, 0) = -4= ^ ^m(±^ r) e me ; m = -1, 0}. 

V 27T 

Thus 

1 /I ? #r 2 \ / Br" 

f 0>± (r,e) = N r a G(- + aT^ ,1 + a,— J exp 7 



■ 



where iV_i and iVp are normalisation constants making the basis orthonormal. 
We shall need the explicit values of iV_i and N . Using the relation 

W V)T {z) = z T+ h- z/2 G Q - v + t, 2r + 1, z^j 

where W is the Whittaker function we get 

POO 

N m ~ 2 = I \g 2 m (±t;r)\ 2 rdr 

x- x W e , a {x) W s>a {x) dx 



w 2 \ \m+a\+l 



2 \B 

where 



oc 



1 ( 1 \ 1 , 

Q = - y-m -a + —J, a= - \m + a\. 
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Combining the identities (T(| 2.19.24.6] 



7T 

x W Btff (x) W„ tCr (x)dx 



sin(27T(j) 

1 (\ 3 

r(i-(x-/i) r(f + v 2 2 

1 (\ 3 

+ VF73 V 2jFi U ~ ° ~ /x> 1; 9 ~ CT ~~ ^ ; 1 

T\± + o-\l)T\\-o-q) V 2 2 



and 



2 F 1 (a,b;c;z) = 

T(c) T(c — a — b) „ . , , . 

F \ F , 7T 2^1 (a, 6; a + b - c + 1; 1 - z) 

I [c — a)l [c — b) 

+ I ^r^P l (! - ^r a ~V\(c - a, c - b;c-a-b + l;l-z) 

r(a)r(6) 

we arrive at the relation 



x 1 W Qia (x)W tl , CT {x)dx 



71 



o ' ' sin(27ro")(/i — g) 

1 1 



X 



r (| - fi - a) T (| - g + a) T (± - p + a) T (§ - g - a) 

Finally we get 



-1/2 



2j V 2n \ r(-i + a+^)r(i-^) 

1/2 



2 ° /sml7ra, . , 
iV n = — \ — - — - Im 



2J \ 2n y r(l + ^)r(l + a-^) 
Let us have a look at the asymptotic behaviour at the origin of the basis 
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functions in the deficiency subspaces M± l . By (|4|) and (0) we have 



g\{±i- r) = a_ li± r~ 1+a + 6_ ll± r 1 ^ + 0(r 1+Q ), 
^ 2 (±2; r) = a ,± r" a + 6 ,± r a + 0(r 2 " a ), 



where 



r(i-a) fB\- 1+a r(-l + a) 



Via) fB\- a , T(-a) 



The coefficients a TO> ±, 6 mi ± are related to the normalisation constants iV m for 
it holds true that 

det M_x = -— (iV_!)- 2 , det M = -- (iV )- 2 . (11) 
1 — a a 

where 

M m = 



Particularly, we shall need the fact that the matrices M_i and M are regular. 

Let us now describe the closure of the operator L. In virtue of the de- 
composition (H) we have 

— v — ^® — 

l = y £ m 

where L m = (L*)^. As it is well known, ^ G Dom(L*) belongs to Dom(L) 
if and only if (tp,L*{p) = (L*ip,(p) for all cp G Af % + Af- t . Thus (L*) m = L m 
for m ^ {-1,0}, and if m G {-1,0} then ip(r) e md G Dom((L*) m ) belongs 
to Dom(L m ) if and only if 



lim rW((p(r),g m (±t,r)) = 

r— >0+ 
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where W(f, g) = {d r f)g — f d r g is the Wronskian. Using the asymptotic be- 
haviour (|I0D and the regularity of the matrix M m we arrive at two conditions 

lim r ^ 0+ (-|m + «| H m+a l^(r) - r^ m+a ^ l d r ip(r)) = 0, 
lim r ^ 0+ (|™ + a\ rl m+Q V(r) - rl m+Q l +1 <9 r ^(r)) = 0, 

which can be rewritten in the equivalent form, 

lim r- 2|m+a|+1 <9 r (r |m+Q| v?(r)) = 0, lim r |m+Q ^(r) = 0. 

But since 

r -|m+ah ( r )| < 1 sup \ x -^\ m +^d x {x\ m+a ^{x))\ 

2\m + a\ (B6 ]o ir [ 

we finally get a sufficient and necessary condition for <p(r) e tmd G Dom((L*) m ) 
to belong to Dom(L), namely 

lim r +a (p(r) = and lim r a (p'(r) =0 if m — —1, 

(12) 

lim r a (f(r) = and lim r a+l ip\r) = if m = 0. 

r^0+ ' r— ►()+ 

This shows that if if) e Dom(L*) = Dom(Z) + M l + M- t then 

^(r,0) = ($}(^)r- 1+a + $^(^)r 1 - Q )e- je + ^(^)r- a + $2(^)r a 
+ a regular part. 
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Let us formally introduce the functionals <3>^ on Dom(L*), 



$1 (V>) 



lim r 

r->0 + 




*2 M 



lim r 




(V>) 



lim r° — 

r-»0 + 2-K 



/ ^(r,0)d0, 



lim r 

r->0+ 




Notice that the upper index refers to the sector of angular momentum while 
the lower index refers to the order of the singularity. If ip & Dom(L) then 
according to (|12[) it actually holds $j(if>) = for j = 1, 2, k = -1, 0. On the 
other hand, if ip G Af l +J\f- l and $*(V0 = for all indices j = 1, 2, = -1, 0, 
then ij; — (this is again guaranteed by the regularity of the matrices M_ x 
and M ). 

Let us introduce some more notation. It is convenient to arrange the 
functionals $^ into column vectors as follows, 



Further, applying the functionals to the basis functions in J\f l -\-H- l we obtain 
four 2x2 diagonal matrices. More precisely, set 



Now it is straightforward to give a formal definition of a self-adjoint ex- 
tension H u of the symmetric operator L determined by a unitary operator 
U : jVj — > M-i. We identify U with a unitary 2x2 matrix via the choice 




(*i,±) 



Then 
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of the orthonormal bases {/-i,±, /o,±} in A/± ? .The self-adjoint operator H u 
is unambiguously defined by the condition: H u C L* and ip G Dom(L*) 
belongs to Dom(H u ) if and only if 

^)eR*n(^ + ^ U J. (13) 



However condition (|T3|) is rather inconvenient and we shall replace it in the 
next section by another one which is more suitable for practical purposes. 

5 Boundary conditions 



To turn (|13| ) into a convenient requirement which would involve boundary 



conditions we shall need the following proposition. Set 



D 



I- a 
a 



There is a one-to-one correspondence between unitary matrices U G U{2) 
and couples of matrices Xi,X 2 G Mat (2, C) obeying 

rank ( ^ \ = 2 (14) 

and 

X*DX 2 = X*DX X (15) 

modulo the right action of the group of regular matrices GL(2, C). The one- 
to-one correspondence is given by the equality 

Ran ( 1 j G Ran ( X ' + 

\X 2 J \ $ 2i+ + $ 2t _U 

Let us note that the equivalence class of a couple (X t , X 2 ) modulo GL(2, 
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corresponds to a two-dimensional subspace in C 4 and hence to a point in the 
Grassmann manifold (^(C 4 ). The complex dimension of <Gi2(C 4 ) equals 4, 
i.e. dini]R <Gj2(C 4 ) = 8. The points of G2(C 4 ) obeying the ("real") condition 
( |j~5| ) form a real 4-dimensional submanifold which is diffeomorphic, according 
to the proposition, to the unitary group U(2). 

To verify the proposition we first show that to any couple (X l5 X 2 ) with 
the properties (ff^), ( fl5| ) there are related unique Y G GL{2, C) and U G U{2) 
such that 



(16) 




where we have set 




/ ALia_i i+ iV-ja-i, 
N a , 



,+ 







V 



iV 6o, + 
Using (|TT|) one easily finds that J is regular and 







\ 

iVoao,- 


N b ,- J 



J 1 



D 
D 




Let us introduce another couple of matrices, V+,V_ G Mat (2, C), by the 
relation 





thus V± = =fz/}($2,±Xi - $i i± X 2 ). It follows that 



V±V± = ( X* X* 




15 



and, consequently, 

v *v_ _ v+ * v+ = ( x * x * W ^ -» D H ^ j = ^ - x;dx 2 

for and .D commute (all of them are diagonal), $* ± = and 

(cf. (|TlD). Owing to the property (|i~5|) we have 

V*V_ = V*V + (17) 
which jointly with the property ([14]) implies that 

Ker V- = Ker V+ = Ker ( V " | = Ker ( Xl | = 0. 

\V + ) \x 2 ) 

The only possible choice of the matrices Y and U satisfying (f0|) is 

Y = Vj\ U = V+VJ 1 . 



The matrix U is actually unitary because of fll?!). 

Conversely, we have to show that any couple of matrices X\, X 2 related 
to a unitary matrix U according to the rule 

X 2 ) \U 

obeys ([i4|) and ([i~5|) . Condition (|i~4]) is obvious since J is regular and condi- 



tion ( |15|) is again a matter of a direct computation. In more detail, since it 
holds 



x*dx 2 - x;dx 1 = ( J [/* ) J* f 



-DO \ U 
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it suffices to verify that 

:)'-•(-)■ 

This concludes the proof of the above proposition. 

Using this correspondence one can relate to a couple Xi,X 2 G Mat(2,C) 
obeying (|T4]) and (|15|) a self-adjoint extension if determined by the condition 

i) G Dom(H) [ ) G Ran | 1 ) . (18) 

V ^ J V x 2 y 

Two couples (X 1; X 2 ) and (X^X^) determine the same self-adjoint exten- 
sion if and only if there exists a regular matrix Y such that (X[,X 2 ) = 
(X{Y, X 2 Y). Moreover, all the self-adjoint extensions can be obtained in this 
way. 

We shall restrict ourselves to an open dense subset in the space of all 
self-adjoint extensions by requiring the matrix X 2 to be regular. In that case 
we can set directly X 2 — I and rename X\ — A. Thus A is a 2 x 2 complex 
matrix satisfying 

DA = A*D. (19) 

The corresponding self-adjoint extension will be denoted H A . The condition 
(lilf ) simplifies in an obvious way. We conclude that H A C L* and ip G 
Dom(L*) belongs to Dom(ii A ) if and only if 

$i(V) =A$ 2 (V), (20) 

and this is in fact the sought boundary condition. 

Matrices A obeying ( |T9| ) can be parametrised by four real parameters (or 
two real and one complex). We choose the parameterisation 

A = , u,v G R, w G C. 

V (1 — a)w v J 
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The relation between A and U reads 



A = + $ lt _U)($ 2t+ + %2,-U)- 1 (21) 

(provided the RHS makes sense). 

The "most regular" among the boundary conditions is $i(V0 = 0, i.e. 
the one determined by A = 0, and the corresponding self-adjoint extension 
is nothing but the standard Aharonov-Bohm Hamiltonian H AB discussed in 
Section |3[ According to (|2lp H AB corresponds to the unitary matrix 



U = -$ x = diag 



6 The spectrum 



r (§ + ^) r(i + « + ^) 



^(2 2b) ^ (2 ^~ 01 2b) 



Let us now proceed to the discussion of spectral properties of the described 
self-adjoint extensions. It is clear from what has been explained up to now 
that everything interesting is happening in the two critical sectors of the 
angular momentum labeled by m = — 1 and m = 0. To state it more formally 
we decompose the Hilbert space into an orthogonal sum of the "stable" and 
"critical" parts, 



H = H s © H c 



where 



H s = l2 ( M +> r dr ) ® C etm(> i n c = L 2 (R + , r dr) ® (C e~ l9 ©CI) 

meZ\{-l,0} 

A self-adjoint extension H A decomposes correspondingly, 

H A = H A \ Hs ®H A \ Hc , 
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and we know that on TC S the operator H A coincides with the standard AB 
Hamiltonian, 

ttA\ _ ttAB\ 

n \h s — ti \ Hs . 

Thus 

a(H A ) = a(H AB \ Hs )Ua(H A \ Hc ) 

and, as explained in Section |3|, 

a(H AB \ Hs ) = {B(2k + 1); k G Z + } U {B(2k + 2a + 1); k G N} 

where the multiplicity of the eigenvalue B(2 k + 1) is infinite while the mul- 
tiplicity of the eigenvalue B{2k + 2a + 1) equals k. On the other hand, 

o(H AB \ Hc ) = {B(2k + 1); k G Z+} U {B(2k + 2a + 1); ke Z+} 

where all the eigenvalues are simple (the first set is a contribution of the 
sector m — — 1 while the second one comes from the sector m = 0). Since the 
deficiency indices are finite the Krein's formula jointly with Weyl Theorem 
||Tl , Theorem XIII. 14] tells us that the essential spectrum a ess (H A \ft c ) is 
empty for any A. Thus the spectrum of H A \fi a is formed by eigenvalues which 
are at most finitely degenerated and have no finite accumulation points. 

Let us derive the equation on eigenvalues for the restriction H A \n c . Let 
A 6 1. In each of the sectors m = —1,0 there exists exactly one (up to a 
multiplicative constant) solution of the equation (L*) m f = \f which is L 2 - 
integrable at infinity (with respect to the measure r dr) and we may take for it 
the function g^(A; r) e tmd (cf. (|])). For a second linearly independent solution 
one may take g^(X; r) e imd provided f3(m, A) £ -Z + (cf. (|)). If (3(m, A) G 
— Z + then a possible choice of a second linearly independent solution is 



r |m+a| H ( p(m, A), 7 (m) 



Br \ / Br 
exp 
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where 



H(P, j,z) = z 1 ^ F{(3 - 7 + 1, 2 - 7,2) 

(cf. ©). 

Thus A is an eigenvalue of -ff A |w c if and only if there exists a vector 
(/!, v) G C 2 \ {0} such that the function 

Mr,0)=^g 2 _ 1 (X;r)e~ l9 + ug 2 (X;r) 

satisfies the boundary condition fl20|). Using again (|4]) and (|7|) one finds that 



where 



o do / V v 7 V o &o / V ^ 



n r(i-q) /i?y 1+a _ r(-i + a) 

a -l - ^71 TT ,0-1 



r(a) i, - r ^ 



r(§ + «-£) V27 ' u r(|-A)- 

This immediately leads to the desired equation on eigenvalues which takes 
the form det A = where 



o,q I \ u I \0 6 



After the substitution 



* = 7j-^> i.e.A=B(l-2*), 
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we get 



T(l-a)T(a) 2 T(a)T(a - 1) (2 

— I u 



T{z)T(z + a) B T(z + a - l)T(z + a) \B 



ra-o)r(-a) / 2 X J - n 



— V 



T(z) 2 \B 

r(a — i)r(— a) , , .. „, 

+ \ } 1 L ( UV - a (l - a) w 2 ) = 0. 
T{z)T(z + a — 1) 

To simplify somewhat the form of the equation it is convenient to rescale the 
parameters as follows, 

fB\ x - a r(a) fB\ a T(l-a) [B . , , , 



2y r(2-«) ' ' V 2 / r(i + «) ' * V 2 

Finally we arrive at an equation depending on three real parameters £,?7,C> 
namely 

i £ y in-? 

T{z) T(z + a) + T(z + a - 1) T(z + a) + T(z) 2 + T(z) T(z + a - 1) ~ 

(23) 

There is no chance to solve equation fl23[) explicitly apart of some partic- 
ular cases. One of them, of course, corresponds to the standard AB Hamil- 
tonian. This case is determined by the values of parameters £ = rj = ( = 
and the roots of fl2"3| ) form the set — Z + U (—a — Z + ). Consider also the 
case when when £ = rj = and ( ^ with the set of roots equal to 
— Z + U(— a — Z + )U{1 — a + (~ 2 }. Comparing the latter case to the former one 
we see that there is one additional root, namely 1 — a + (~ 2 , which escapes 
to infinity when ( — > 0. 

In the last particular case one can also consider the limit ( — > oo. More 
generally, suppose that det A / 0, i.e. £77 — £ 2 7^ 0, replace A with t A in (|20f) 
and take the limit t — > 00. The limiting boundary condition reads 

$ 2 (V0 = 
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and the corresponding self-adjoint extension which we shall call H°° is one 
of those omitted when we restricted ourselves to an open dense subset in 
the space of all self-adjoint extensions (regarded 4-dimensional real 

manifold). Equation (^3|) reduces in this limit to the equation 

' — = (24) 



T{z)T{z + a 
with the set of roots — Z . U (1 — a — Z_ 



Another case when equation fl23[) simplifies though it is not solvable ex- 
plicitly is C — 0. This is easy to understand since if C = then the matrix 
A is diagonal and the two critical sectors of angular momentum do not in- 
terfere. This is reflected in the fact that the equation ( p3|) splits into two 
independent equations, 

If 1 77 

+ " 7T = 0, — r + —At = 0. 



T(z) T(z + a-l) T(z + a) T(z) 

Let us shortly discuss the dependence of roots of equation ( |2~3|) on the 



parameters f, 77, C- Since the derivative of the LHS of fl23|) with respect to z 
and with the values of parameters (f , 77, () = (0, 0, 0) equals 

(—\\ m m\ (— l) m m! 

7^ for z = —m, and — 7^ for z = —m — a, 



r(— m + a) T(— m — a) 

where m G Z + , the standard Implicit Function Theorem (analytic case) is 
sufficient to conclude that the roots are analytic functions in f , 77, £ at least in 
some neighbourhood of the origin (depending in general on the root). Let us 
denote by zi tTn (£,r],Q an d £2,m(£,?7,C) the roots of (f2~3|) regarded as analytic 
functions in f , 77, £ and such that £i, m (0, 0, 0) = — m and 2; 2im (0, 0, 0) = —a — 
m, with m G Z + . A straightforward computation results in the following 
power series truncated at degree 4. 
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Set 



h° m (z) 


= E 










h l m {z) 


6 




= -2 



J 



J'=l J 



J'=l J 

where 7 is the Euler constant, ^(z) = r'(;z)/r(,z) is the digamma function 
and £ is the zeta function. Then 

ZlM) = _ m+ (-'>"" ; {+ 7 g-'-^;) . 



+ 



m! T(— 1 — m + a) (m!) 2 T(— 1 — m + a) 2 
(-l) m+1 (3 - m + a) 2 + hji-l -m + a)) 

2(m!) 3 r(-l -m + a) 3 



| (-IT (1 + m-a) 2 
m! r(— 1 — m + a) 

+ eWrFT^W (4 - m + a) (25) 

x (4 ^(-1 - m + a) 2 + 3 ^(-1 - m + a)) 

+ / i 2 t (-l-m + a))e 4 

3 - 2 (1 + m - a) /^(-m + a) 2 2 

(m!) 2 r(-l-m + a) 2 ^ ' 
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m\T{—m — a) (m!) 2 T(— m — a) 2 

, (-ir 1 (3^(- m -a) 2 Ui(-m- a )) T?3 

2 (m!) 3 T (— m — a) 3 

I (-^ + 2 
m! r(— m — a) 

+ R , ,u r / u(4CU-«) (26) 

x (4/i^(— m - a) 2 + 3 h 2 m (-m - a)) 
+ h 2 m (-m - a)) r/ 4 
, l-2(m + l)^(-m-a) 2 2 + 
(m!) 2 r(— m — a) 2 

A similar analysis can be carried out to get the asymptotic behaviour of 
roots for £, rj, ( large. To this end assume that £77 — £ 2 7^ and set 

£ v H C 



Notice that £V — (' 2 = (£77 — C 2 ) -1 - Equation ( f23| ) becomes 



T(z)T(z + a) V{z + a- l)T(z + a) T(z) 2 T(z) T(z + a - 1) 

(27) 



Roots of flgjj ) are analytic functions in 7/, at least in some neighbourhood 
of the origin. Again, it would be possible to compute the beginning of the 
corresponding power series and to derive formulae similar to those of (|2ol) , 
(|26|) but we avoid doing it here explicitly. 

Instead we prefer to plot two graphs in order to give a reader some im- 
pression about how the eigenvalues may depend on the parameters, i.e. on 
the boundary conditions. In each graph we choose a line in the parameter 
space, {(££, rjt, (t) £ R 3 ; t £ M}, and we depict the dependence on t of several 
first eigenvalues for the corresponding self-adjoint extension restricted to Tt c 
(see (|2"2"D for the substitution). In the both graphs we have set a = 0.3 and 
B = l. 
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Probably the most complete general information which is available about 
solutions of equation fl2"3| ) might be a localisation of roots of this equation 
with respect to a suitable splitting of the real line into intervals. Let us 
choose the splitting into intervals with boundary points coinciding with the 
roots of equation (^J|). To get the localisation let us rewrite equation (p3|) , 
equivalently provided z ^ — Z + U (1 — a — Z + ), as follows 



Y{z - 1 + a) A ( T(z) , 2 



T{z) J \T(z + a) 

Put 

T(z-l + a) 



K hr^ + lh^. (28) 



FJz) 



T(z) 



so that equation fl28|) can be rewritten as 

(F a (z) + {F x . a {z + a) + V ) = C 2 . (29) 

It is easy to carry out some basic analysis of the function F a (z). We 
have F a '(z) = F a (z) (i/j(z — 1 + a) — i/i(z)). One observes that F a (z) > 
for z G ]1 — a, +oo[ U ([J meZ+ ] — a — m, — m[ j, and F a (z) < for z G 
|J meZ+ ] — m, 1 — a — m{ , and in any case F a '(z) < 0. In the former case this 
follows from the fact that ip( z ) is strictly increasing on each of the intervals 
]0, +oo[ and ] — m— 1, — m[, with m G Z + . In the latter case this is a 
consequence of the identity 



7r sin 7ra 



oo g-(l-z)* ^ _ 



^(2 — 1 + at) — 1p(z) — -r-, r — r-r~/ rr + / : —, dt. 

sm{irz) sin(7T(2 + a)) J 1 — e 1 

Moreover, 

lim F a (z) = 0, lim F a (z) = ±00 and F a (— m) = for m G Z + . 
This also implies that Fi^ a (z+a) > for z G ]0, +oo[ U ^UmfEZ+ ] — 1 — m > ~ a 
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and Fi_„0) < for z G \J m& 
0, and 



a — m, —m[ , in any case Fi- a '(z + a) 



< 



lim Fi_ a (z + a) = 0, 
and Fi_ Q (— a — m) = for m G Z_| 



lim Fi_ a (z + a) = ±oo, 



With the knowledge of these basic properties of the function F a (z) it 
is a matter of an elementary analysis to determine the number of roots of 
equation (|29|) in each of the intervals ]1 — a, +oo[ , ] — m, 1 — a — m[ and 
] — a — m, —m[ , with m G Z + . The result is summarised in the following 
tables. 



interval ]1 — a, +oo[ 


conditions 


number of roots 




v > o 


c 2 >^ 


1 


e>o 


?7 > 


C 2 <£r/ 







-r(l - a) < 77 < 


no condition 


1 


£>o 


r? < -r(i - a) 


no condition 





£<o 


v > o 


no condition 


1 


e<o 


-r(l - a) < r? < 




1 


e<o 


-r(l - a) < 7] < 




2 


e<o 


77 < -r(l - a) 


c 2 >^ 





£<o 


77 < -r(i - a) 




1 



interval ]0, 1 — a[ 


conditions 


number of roots 




77 > -r(i-a) 





£<o 


77 < -r(i-a) 


1 


e>o 


77 > -r(i - a) 


1 


£>o 


r) < -r(i - a) 


2 
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intervals ] — a — m, —m{ , m G Z + 


conditions 


number of roots 


£>o 


77 < 





e>o 


T] > 


1 


£<o 


r? < 


1 




r/ > 


2 




intervals ] — 1 — m, —a — m[ , m G Z + 


conditions 


number of roots 




r? > 







1] < 


1 




77 > 


1 




7/ < 


2 



This is to be completed with the simple observation that 1 — a is a root of 
( |23|) if and only if r\ = —T(l — a), and —to, with m G Z + , is a root if and 
only if £ = 0, and finally — a — m, with m 6 Z + , is a root if and only if r] = 0. 

Let us note that this localisation is in agreement with a general result 
according to which if A and B are two self-adjoint extensions of the same 
symmetric operator with finite deficiency indices (d, d) then any interval J C 
]R not intersecting the spectrum of A contains at most d eigenvalues of the 
operator B (including multiplicities) and no other part of the spectrum of B 



[|I2|, §8.3]. Thus in our example if J is an open interval whose boundary points 
are either two subsequent eigenvalues of H°° or the lowest eigenvalue of H°° 
and —00 then any self-adjoint extension H A has at most two eigenvalues in 
J. 



7 Concluding remarks 

The above discussion does not exhaust all questions related to the system 
under consideration. One may ask, for instance, how the state of such a 
particle evolves under an adiabatic change of parameters. In particular, 
since the model exhibits eigenvalue crossings, one may expect that there 
are parameter loops exhibiting a nontrivial Berry phase. Another question 



27 



concerns the physical meaning of our idealized model. More specifically, one 
is interested in which sense the model Hamiltonian can be approximated by 
those with smeared flux and a regular interaction. We leave these problems 
to a future publication. 
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Figure 2: The Hamiltonian is determined by the boundary conditions corre- 
sponding to the parameters C) — (0.95 1, —0.25 1, 0), a = 0.3, B — 1. 
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